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of  danse,  partially  Ionised,  reacting  plasmas  has  bean  carried  ant.  This 
rapoitwjmartm  the  work  *1ch  Is  fully  described  in  14  journal  articles, 
and  reports.  \ 
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A  RESEARCH  PROGRAM  FOR  THE  THEORETICAL  STUOY 
OF  THE  THERMODYNAMIC  AND  TRANSPORT  PROPERTIES  OF  DENSE  PLASMAS : 

FINAL  REPORT* 

F.  J.  Rogers,  H.  E.  DeWItt,  and  D.  B.  Boercker 
University  of  California,  Lawrence  Livermore  National  Labortory 
Livermore,  California  945SO 

I.  INTRODUCTION 

Over  the  past  five  years,  we  have  pursued  a  program  of  theoretical 
research,  partially  funded  by  the  Office  of  Naval  Research,  on  the  Equation 
of  State  (EOS)  and  transport  properties  of  dense,  partially  Ionized,  react* 
Ing  plasmas.  Although  originally  motivated  by  the  failure  of  older,  simpler 
theories^  to  describe  plasmas  produced  In  new.  Innovative  power  sources, ^ 
this  program  developed  Into  a  far-ranging,  multi-faceted  study  of  the  basic 
properties  of  such  complicated  systems.  Our  research  efforts,  nevertheless, 
can  be  divided  Into  three  basic  categories:  1)  Theory  of  Equation  of  State 
and  Static  Structure;  2)  Electrical  Conductivity  and  Dynamic  Properties,  and 
3)  Computer  Simulation  Studies.  While  these  categories  are  fairly  well  de¬ 
fined,  they  are  very  closely  related,  and  results  from  each  Impact  on  results 
from  the  others.  Since  all  of  the  most  significant  results  have  already 
appeared  In  published  form,  the  progress  made  In  each  of  these  areas  Is 
described  only  briefly  In  the  following  paragraphs.  More  detailed  descript¬ 
ions  can  be  found  In  the  references  listed,  by  year  of  publication.  In 
Section  HI  of  this  report. 


*flerinp«f7of5Mf  under  the  auspices  of  the  U.S.  Department  of  Energy  by 
Lawrence  Livermore  National  Laboratory  under  contract  M-7405-Eng-48  and 
supported  In  part  by  the  Office  of  Naval  Research. 


II.  MAJOR  ACCOMPLISHMENTS 


A.  Theory  of  Equation  of  State  and  Static  Structure 

A  significant  part  of  our  work  for  ONR  Involved  the  development  of  a 
fundamental  theory  of  the  Ionization  balance  and  equation  of  state  of  dense, 
reacting  plasmas.  Under  ONR  funding,  we  were  able  to  extend  earlier  work  on 
hydrogen  and  helium  to  Include  heavier  e1ements[81a,bj.  A  detailed  tabular 
EOS  (with  a  polynomial  fit  provided)  was  constructed  for  argon[80a].  Similar, 
but  less  detailed,  calculations  were  also  carried  out  for  hydrogen,  helium, 
and  xenon[83a].  Particular  emphasis  was  placed  on  the  rare  gases,  because 
they  were  anticipated  to  be  the  working  fluids  In  the  ARTEC  MHO  power 
generator/ 

To  facilitate  these  calculations  we  had  to  develop  effective  potentials 
to  describe  the  Interactions  between  free  electrons  and  Ions  with  up  to 
fifty-five  bound  electrons.  The  model  we  constructed  Is  mathematically 
slwple,  but  It  still  retains  the  shell  structure  of  the  1ons[81c]. 

Another  Important  feature  of  dense  plasmas  that  we  studied  was  their 
equilibrium  static  structure.  This  Is  most  easily  described  In  terms  of  the 
static  structure  (or  font)  factors  of  the  plasma.  These  functions  are  related 
to  the  Fourier  transform  of  the  pair  distribution  function,  which  gives  the 
probability  of  finding  a  pair  of  plasma  particles  at  a  given  separation.  In 
weakly  coupled  plasmas  the  pair  correlations  are  well  described  by  the 
Debye -Hucke!  theory/  but  for  dense  plasmas,  other  methods  must  be  used. 

Maong  the  most  successful  of  these  Is  an  Integral  equation  approximation, 
known  as  the  Ifrpernetted  Cheln(MC)  Equation/  This  equation  has  been  used 
to  described  the  classical  one-component  plasma  end  molten  salts  at  moderate 
coup11ng[80b,c].  For  very  strong  coupling,  however,  the  NNC  equation  must  be 
improved.  Me  have  succeeded  in  doing  this  by  blending  the  NNC  equation  with 
another  Integral  equation,  called  the  Fercus-Tevlck  equation.  In  a 
thermodynamically  consistent  way.  The  blended  Integral  aquation  gives  the 
accuracy  of  very  expensive  computer  simulations  and  Is  no  mare  difficult  to 
solve  than  the  NNC  oquet1on(Mb,c]. 
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For  real  plasmas  consisting  of  electrons  Interacting  with  Ions  through  an 
attractive  Coulomb  potential,  direct  application  of  these  classical  Integral 
equation  techniques  Is  Impossible,  because  quantum  mechanics  Is  required  to 
prevent  the  collapse  of  the  system.  However,  If  the  density  Is  not  too  great, 
the  essential  quantum  mechanics  can  be  Incorporated  Into  a  pseudopotentlal 
which  Is  finite  at  Infinitesimal  distances,  but  behaves  like  the  Coulomb 
potential  at  large  distances.  Using  results  from  our  equation  of  state 
theory,  we  have  found  a  suitable  method  for  determining  these  numerical 
pseudopotentlal s  from  our  analytic  effective  potentials  described  above. 

Having  accomplished  this,  we  can  use  the  pseudopotentlal s  In  the  HNC  equation 
to  obtain  both  the  static  structure  and  the  equation  of  state  properties  of 
real  plasmas  from  the  same  theory[64a,d]. 

B.  Electrical  Conductivity  and  Dynamic  Properties 

Our  first  major  success  In  this  area  was  the  calculation  of  the 
electrical  conductivities  of  argon  and  xenon  plasmas  using  a  Chapman-Enskog 
solution6  to  the  Bould,  Ml 11 lams,  DeWltt(GMD)  kinetic  equation.7  These 
calculations  corrected  two  of  the  Important  shortcomings  of  the  usual 
Spltzer-type  theories.  First,  the  relevant  statically  screened  cross  sections 
were  calculated  from  their  correct  quantum  expression.  Second,  scattering 
from  the  electronic  shells  of  the  Ions  and  from  neutrals  was  Included  through 
the  use  of  analytic  effective  potentials  developed  as  part  of  this 
ONR-sponsored  program.  These  two  corrections  greatly  Improved  the  agreement 
with  Soviet  experiments  In  moderately-coupled,  partially  Ionized 
plasmas(81d].  Dynamic  screening  effects  were  also  studied  In  the  context  of 
the  QUO  kinetic  equat1on[81e],  but  the  results  for  systems  with  Important 
short-range  contributions  to  the  potential  were  not  as  accurate  as  the  Zlman 
formula8  or  the  simple,  statically  screened  resu1ts£82a]. 

Ha  have  also  studied  the  relationship  between  the  Zlman8  and  AMD 
theories  of  conductivity,  and  we  shewed  that  those  apparently  different 
approaches  can  be  obtained  from  a  single,  unified  approach  based  on  the 
8reen-Kube  theory8  for  the  electrical  conductivity.  The  net  result  of  this 
new  approach  Is  a  simple  expression  for  tbs  electron-ion  collision  frequency 
■bleb  depends  only  on  the  static  structure  factors  for  the  p1emae£«b}.  Per 


modest  densities,  this  theory  gives  Improved  agree— nt  over  the  theory  of 
Spitzer  for  aval lable  experl— nts  on  hydrogen  plasms  and  It  also  gives  good 
agree— nt  with  Molecular  dynamics  simulations  of  strongly-coupled  plasmas[82c]. 

The  new  theory  described  above  had  the  additional  advantage  that  It  could 
be  extended  to  non-zero  frequencies  to  study  the  conductivity  as  a  function  of 
frequency.  Once  again,  comparison  with  — lecular  dynamics  results  for 
frequencies  below  about  half  the  plasma  frequency  gave  excellent  agreement. 
Extending  the  calculations  beyond  the  plasm  frequency,  predicted  a 
substantial  reduction  of  the  plasm  — Isslvlty  by  short-range 
corre1at1ons[83b]. 

Our  Interest  In  dynamic  conductivities  required  a  study  of  the  more 
general  dynamic  strucure  properties  of  plasms.  A  comparison  of  several  model 
kinetic  equations  for  the  dynamic  structure  factors  has  shown  that  the  theory 
Llnnebur  and  Duderstadt10  accurately  describes  the  electron  peak  In  the 
charge-charge  structure  factor  generated  In  simulation  stud 1es[ 83c] 

C.  Computer  SI— latlon  Studies 

He  have  conducted  very  detailed  and  extensive  nu— rlcal  si— latlon 
studies  of  the  thermodynamic  properties  of  strongly  coupled  plasms.  In  this 
work,  the  Monte  Carlo  — thod11  was  used  to  compute  statistical  averages  over 
hundreds  of  thousands  of  configurations  for  several  hundred  Ions  in  themo- 
dyn— 1c  equilibrium.  The  problems  —  studied  Included:  1)  the  classical  one 
component  plasm  (OCP),  consisting  of  point  ions  In  a  uniform  background; 

2)  Ionic  mixtures;  3)  point  Ions  In  a  polarizable  background,  and  4)  the 
evaluation  of  the  Ion  mlcrefleld  distribution  for  each  of  the  previous  three 
system. 

for  the  OCf.  —  were  able  to  obtain  the  beet  estl— te,  to  date,  of 
liquid-solid  phase  trans1t1oe(82d],  and  —  obtained  very  accurate  results  for 
the  pair  distribution  function  and  structure  factors.  Those  latter  results 
were  very  Instrumental  In  our  study  of  the  blended  Integra!  equation 
previously  described,  free  ear  wet*  —  Ionic  mixtures  m  established  a  staple 
and  accurate  nixing  rule  for  date-doing  the  MS  of  a  mixture  free  the 
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properties  of  the  Indvldual  component s[84e].  In  order  to  model  more  realistic 
systems,  we  used  linear  response  theory  to  model  the  electrons  as  a 
polarizable  background  which  screened  the  Ion-ion  Interactions,  and  applied 
this  model  both  to  one-component  systems  and  to  Ionic  m1xtures[83d]. 

The  observation  of  spectral  lines  emitted  from  plasmas  yields  Important 

Information  about  the  sate  of  the  plasma.  As  a  consequence,  we  modified  our 

simulation  codes  to  calculate  the  central  Ingredient  to  most  theories  of 

plasma  line  broadening,  the  distribution  of  electric  microfields  around  a 

radiating  ion.  Our  calculations  provided  the  benchmark  checks  for  a  new  and 
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very  efficient  theory  for  the  microfield  called  APEX.  This  new  theory 
accurately  describes  the  microfield  distributions  In  all  three  of  the  plasma 
models  described  above[84f]. 

III.  References  to  ONR  Sponsored  Research 

In  this  section  we  list  the  work  we  published  under  ONR  sponsorship  from 
1980  to  1984  by  year  of  publication.  These  Include  14  articles  In  refereed 
journals,  4  Internal  reports,  and  2  Invited  papers. 

1980 

a.  "Equation  of  State  for  Self-Excited  MHO  Generator  Studies",  F.  J.  Rogers, 
N.  Ross,  6.  L.  Haggle,  and  l.  K.  Hong,  LINL  Report,  fJCI0-18SS7.  Feb., 
1980. 

b.  "Improved  Equation  of  State  for  the  Classical  One-Component  Plasma".  M. 

L.  Slattery,  6.  0.  Ooolen,  and  H.  E.  OeHItt,  Phys.  Rev.  A  21,  2067  (1980). 

c.  "A  HNC  Study  of  Asymmetrically  Charged  Hard  Spheres",  F.  J.  Rogers, 

J.  Chen.  Phys.  73,  8272  (1980). 

i&l 

a.  "Equation  of  State  of  Donee,  Partially  Degenerate,  Reacting  Plasmas", 

F.  J.  Regers,  Phys.  Rev.  A  24,  1531  (1981). 

b.  "Comparison  of  Two  Equatlon-of -State  Models  for  Partially  Ionized 
Aluminum i  Zel'dovlch  and  Raiser's  Model  Versus  the  Activity  Expansion 
Code",  R.  J.  Merrach  and  F.  J.  Rogers,  j.  Appl.  Phys.  82,  8882  (1981). 

e.  "Analytic  ftoctren-Ion  effective  Potentials  for  2  <  88",  F.  3.  Rogers, 
Phys.  Rev.  A  23,  1808  (1981). 
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d.  "Electrical  Conductivity  of  Dense  Plasmas",  F.  J.  Rogers,  H.  E.  DeWItt, 
and  0.  B.  Boercker,  Phys.  Lett.  82A,  331  (1981). 

e.  "Dynamic  Screening  Effects  on  Electrical  Conductivity  In  Dense  Plasmas", 
D.  B.  Boercker,  F.  J.  Rogers,  and  H.  E.  DeWItt,  LLNL  Report, 
UCID-18574-81-2,  Sep.,  1981. 

1982 

a.  "Dynamic  Screening  of  Electron-Ion  Collisions",  F.  J.  Rogers  and  D.  6. 
Boercker,  Bull.  Amer.  Phys.  Soc.  27,  1105,  (1982). 

b.  "Electron  Collision  Frequency  In  Plasmas",  D.  B.  Boercker,  F.  J.  Rogers, 
and  H.  E.  DeWItt,  Phys.  Rev.  A  25,  1623  (1982). 

c.  "Electrical  Conductivity  of  Hydrogen  Plasmas",  D.  B.  Boercker,  LLNL 
Report,  UCID-18574-81-4,  Mar.,  1982. 

d.  "N-Dependence  In  the  Classical  One-Component  Plasma  Monte  Carlo 
Calculation",  W.  L.  Slattery,  G.  D.  Ooolen,  and  H.  E.  DeWItt,  Phys.  Rev. 

A  26,  2255  (1982). 

1983 

a.  "Ionization  Equilibrium  and  Equation  of  State  In  the  Solar  Interior", 

F.  J.  Rogers,  Proceedings  of  the  Conference  on  Solar  Seismology  From 
Space,  Aug  17-18,  1983. 

b.  "Frequency-Oependent  Electrical  Conduclvlty  of  Strongly-Coupled  Plasmas", 
D.  B.  Boercker,  LLNL  Report,  UCI0-18574-83-2,  Dec.,  1983. 

c.  "Dynamic  Structure  Factors  In  Two-Component  Plasmas",  R.  Cauble  and  0.  8. 
Boercker,  Phys.  Rev.  A  28,  944  (1983). 

d.  "Some  Approximate  Microfield  Distributions  for  Highly  Ionized  Plasmas:  A 
Critique",  C.  A.  Igleslas,  C.  F.  Hooper,  Jr.,  and  H.  E.  DeWItt,  Phys. 

Rev.  A  28,  361  (1983). 

im 

a.  "Integral  Equation  Method  for  Partially-Ionized  Plasmas",  F.  J.  Rogers, 
Phys.  Rev.  A  29,  868  (1964). 

b.  "New,  Thermodynamically  Consistent  Integral  Equation  for  Simple  Fluids", 
F.  J.  Regers  end  D.  A.  Young,  Phys.  Rev.  A  30,  999  (1984). 

c.  "A  New,  Accurate  Statistical  Mechanical  Theory  for  Simple  Fluids",  F.  J. 
Regers  and  D.  A.  Young,  Phys.  Lett.  102A,  303  (1984). 

d.  "Statistical  Mechanics  of  Dense  Plasmas  and  Implications  for  the  Plasma 
Polarization  Shift",  F.  J.  Regers,  In  Proceedings  of  the  Seventh 
International  Conference  on  Spectral  Line  Shapes,  Aussols,  France,  1984 
(In  press). 
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e.  "Statistical  Mechanics  of  Light  Elements  at  High  Pressure  VII:  A 
Perturbative  Free  Energy  for  Arbitrary  Mixtures  of  H  and  He",  W.  B. 
Hubbard  and  H.  E.  DeWItt  (to  appear  in  Ap.  J. ) . 

f.  "Low-Frequency  Electric  Microfield  Distributions  In  Plasmas",  C.  A, 
Igleslas,  H.  E.  DeWItt,  J.  L.  Lebowftz,  and  D.  MacGowan,  (to  appear  in 
Phys.  Rev.  A). 
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